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1 Introduction 

The study of non-relativistic two-body systems consisting of a heavy quark- 
antiquark {QQ) pair has a long tradition since the discovery of the J /ip and the 
work of Appelquist and PoUtzeru, which had shown that, owing to asymptotic 
freedom, non-relativistic quantum mechanics should apply to heavy QQ sys- 
tems. Thus, heavy QQ systems share many features of the positronium and are 
an ideal tool for studying the interplay of perturbative and non-perturbative 
physics. It is clear that non-perturbative effects become suppressed more ef- 
ficiently for increasing heavy quark mass. To get a more quantitative feeling 
about this, it is helpful to consider the scales that are relevant to the dynam- 
ics of a positronium-like QQ pair with relative velocity w ^ 1, neglecting for 
now the influence of non-perturbative effects. There are three: the mass, mq 
(relevant to processes with high virtuality such as annihilation), mqv (which 
sets the size of the system and the relative momentum), and rnqv^ (which is 
of order of the binding energy) . To assess the importance of non-perturbative 
effects, one can compare these scales with the typical hadronic scale Aqcd- 
The tools needed to describe the QQ pair depend on the relation of Aqcd to 
the three scales just mentioned. 

In this talk I consider heavy QQ systems, where mg S> rnqv 3> itiqv^ 3> 
Aqcd, and I call such systems "perturbative'H. In fact, only a top-antitop {tt) 
pair can clearly satisfy this condition, and it satisfies it even for arbitrarily 
small velocities because of the large top decay rate. (I will come back to this 
point in Sec. ^) For a bottom-antibottom (bb) pair, the condition is, if at all, 
only valid for the bottomonium ground state, whereas for a charm-anticharm 
pair it can certainly never be satisfied. For bb pairs, on the other hand, one 
can construct quantities, such as moments of the total cross section of hadrons 

*This write-up was updated for some relevant publications that have appeared after the time 
of the conference. 
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containing a bottom and and antibottom quark in e^e^ annihilation where the 
condition is fairly well satisfied. Perturbative QQ systems are the ones that 
are best understood. They allow for a rather precise quantitative description 
from first principles QCD, i.e. without a sizeable model-dependence. This 
makes perturbative QQ systems in principle an ideal tool to determine QCD 
parameters such as the heavy quark mass uiq and the strong coupling as- 

In this talk I review the recent achievements in the description of pertur- 
bative non-relativistic QQ systems using the concept of effective field theories. 
These developments have allowed a systematic determination of corrections of 
order with respect to the non-relativistic limit for heavy quark production 
(or annihilation) in e~^e~ annihilation and opened up the access to the de- 
termination of even higher order corrections. In Sec. ^ I will briefly motivate 
why the concept of effective field theories is so helpful, and in fact necessary 
to understand perturbative QQ systems quantitatively, and in Sec. ^ I discuss 
basics of the construction of an effective theory for non-relativistic perturba- 
tive QQ systems ignoring non-perturbative effects altogether. In Sec. ^ the 
calculation of the total QQ production cross section in e^e~ annihilation close 
to threshold at NNLO is sketched, while in Sec. |^ the need to introduce heavy 
quark mass definitions that are adapted to the non-relativistic framework is 
explained. Section ^ contains the application of the NNLO cross section to 
the top mass determination at a future Icpton pair collider. In Sec. ^ I briefly 
mention further applications and give the conclusions. 



2 Why do we need an EfTective Theory? 

The lowest order description of the perturbative QQ pair in the non-relativistic 
expansion is provided by the well-known Schrodinger equation 

Because all relevant scales are larger than Aqcd, the coupling is weak and 
the potential can be calculated perturbatively (similar to QED). Equation (^ 
provides a resummation of the most singular terms oc (as/u)" in the elastic 
QQ QQ scattering amplitude for u ^ in full QCD. There is nothing 
more we have to do in the non-relativistic limit. However, a problem arises 
when we try to calculate relativistic corrections. As an example, let us con- 
sider the Darwin term SHd — |^^(5*^'^)(r). The first order corrections to the 

energy levels can be determined trivially as SEn — {n\H]:)\n) ^ |-0n(O)p, using 
Rayleigh-Schrodinger time-independent perturbation theory, where n stands 
for the appropriate set of quantum numbers. For the description of production 
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and decay rates, however, we also need the corrections to the wave function at 
the origin. Formally the answer reads 



(2) 

Interestingly the expression is UV-divergent because G(0, r, E) contains terms 
proportional to 1/r and Inr for r ^ 0. (The power divergence can already 
be found in the free Green function Gheo{0,r, E) — ^^e^P"", p — yjEmq.) 
The same divergence also exists in the energy levels, but in the next order of 
perturbation theory. 

Obviously, the UV-divergence reflects the fact that the non-relativistic 
Schrodinger equation is only an approximation, which is valid for small mo- 
menta of order rnqv and smaller. This means that the summation over the 
intermediate states with very large energies of order mg and larger is incorrect. 
The traditioual approach to avoid this problem has been promoted by Bethe 
and SalpeterS and simply avoids the non-relativistic approximation completely. 
This is achieved by formulating the problem in terms of fully relativistic func- 
tional equations of Green functions. The Bethe-Salpeter approach, however, 
has two disadvantages: first, the corresponding equations are quite difficult 
to solve (because it is up to the person who solves them to find out what 
exactly has to be calculated at a specific order in the non-relativistic expan- 
sion) and, second, there is no systematic way to separate perturbative and 
non-perturbative effects. 



G(0,0,£;) 



En — E 



V^n(O) . 



3 The EfTective Field Theory Approach 

The effective field theory (EFT) approach avoids this dilemma. The EFT 
approach for non-relativistic perturbative QQ pairs had already been initiated 
some time ago by Caswell and Lepage □, but only recently has it been fully 
elaborated. To construct the EFT one has to separate those physical degrees of 
freedom (d.o.f.'s) that can become on-shell for the non-relativistic QQ pair from 
those that can only be off-shell. This separation is, of course, not unambiguous. 
So, the separation between what we count as on-shell and off-shell is dependent 
on our choice of the regularization scheme. The idea is to integrate out the 
off-shell d.o.f.'s, in much the same way as particles with masses that are much 
larger then the available energies can be integrated out in a given process. In 
order to describe production and annihilation processes, it is also necessary to 
carry out the same program for external currents. Once the EFT Lagrangian 
is obtained one can derive equations of motion for Green functions involving 
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the QQ pair. Here, the UV-divergence mentioned before does not arise from 
the fact that the EFT is defined within a regularization scheme. The resulting 
dependences on the regularization parameter are cancelled once the Wilson 
coefficients of interactions and the external currents in the EFT (which contain 
the contributions from the off-shell d.o.f.'s) are taken into account. 

For the case of the non-relativistic QQ pair this program is non-trivial 
because we have to distinguish between on- and off-shell components of the 
single heavy quark d.o.f. in full QCD. The same it true for the gluon^ d.o.f. 
in full QCD. A very useful way to identify the relevant d.o.f.'s for the non- 
relativistic QQ system is to identify the momentum regimes needed to set up 
an asymptotic expansion of Feynman diagrams (in full QCD) that describe 
elastic QQ QQ scattering or QQ production and annihilation for centre- 
of-mass energies very close to 2mQ. It is obvious that, in general, one cannot 
naively carry out the expansion before integration, because the resulting terms 
might not lead to finite results. As a simple mathematical example we might 
consider that we want to determine the expansion of the parameter integral 
/(o) = /g arctan(l — x)/{x^ + a?)dx for a 0. Naive expansion before 
integration will lead to divergent expressions. Here we have to distinguish 
between the regions where a: ~ a and x ^ 1. We might introduce a cutoff A that 
separates the regions such that a ^ A ^ 1, but we could also use an analytic 
regularization scheme. In the region x ^ a we can now conveniently expand 
the arctan for small x, and for a; ~ 1 we can expand the whole integral in a. 
The integrations can now be done rather easily. Adding both contributions 
and expanding properly in A and a, the cutoff dependence cancels and leaves 
us with the correct expansion in a. 

The corresponding momentum regions relevant to a non-relativistic QQ 
pair are more involved, because there are three widely separated scales one 
has to take into account: toq, moij and raqv^ . The regions were first fully 
identified by Beneke and SmirnovB: (k^,k) ~ (mQ.mQ) ("hard"), (fc^jfc) ~ 
{mQV,mQv) ("soft"), {k",k) ^ {niQv'^^mQv) ("potential") and (fc°,fc) ~ 
{niQV^, niQV^) ( "ultrasoft" ) . The energy and momentum components can scale 
differently with v because the choice of a cm. system for the QQ pair breaks 
Lorentz invariance. All these regions can be populated by quark and gluonic 
d.o.f.'s. For the QQ system, where the available energy is of order mQV^ , 
only quarks in the potential regime have a chance to be on-shell. For gluons 
this can only happen in the ultrasoft regime. In order to get to the EFT we 
have to integrate out all the other d.o.f.'s. There have been different attempts 
and philosophies to achieve this. Pineda and SotoQ have proposed a two-step 
procedure, where first the d.o.f.'s in the hard regime are integrated out at 

t We assume that all light quarks are massless, and we count them as gluons. 
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the scale_mQ, leading to "non-relativistic QCD" (NRQCD), the theory formu- 
lated inlj. In a second step, the other ofF-shell d.o.f.'s are integrated out at 
a scale of order rnqv, leading to "potential NRQCD" (pNRQCD). This two- 
step procedure has the feature that it contains two separate matching scales. 
Manohar and collaborators QIj propose a single-step procedure where the ofF- 
shell d.o.f.'s are integrated out at the scale mg, except for soft gluons. The 
resulting EFT is called "velocity NRQCD" (vNRQCD). The crucial issue of 
vNRQCD is that it is regulated by a regularization velocity ly in dimensional 
regularization. This regularization allows the scaling of the boundaries of the 
different momentum regions with v at the same time. A common feature of 
both EFTs is that their Lagrangians contain four-quark {QQ){QQ) interac- 
tions that depend non-locally on the quark three momenta, but locally on the 
quark energies. These interactions represent a generalization of what is called 
"potential" in the Schrodinger equation (|l|). We note that, in the EFT for 
perturbative QQ pairs, these potentials are purely short-distance objects that 
can be calculated perturbatively and do, in contrast to potential models, not 
contain any confining contributions. 

Setting aside the problem of properly resumming logarithms of the ve- 
locity □ (which is acceptable because the corresponding logarithms are not 
dominant numerically for the practical applications involving ti and bb pairs) 
both formulations of the EFT agree on how to calculate the relativistic correc- 
tions to QQ production. From now on, I will therefore only talk about "the 
EFT for perturbative non-relativistic QQ pairs" . Based on the velocity scaling 
of the different momentum regimes, it is possible to unambiguously determine 
to which power in v an individual interaction in the Lagrangian of the EFT 
scales. This "power counting" allows an unambiguous identification of the in- 
teractions (and the contributions in their Wilson coefficients) that have to be 
taken into account to describe a process involving the non-relativistic QQ pair 
at a certain order in v (i.e. in the non-relativistic expansion). One well-known 
fact that can be derived from formal power counting is that the Coulomb inter- 
action scales with the same power in v as the kinetic energy, i.e. the Coulomb 
interaction cannot be treated as a perturbation. As a particularly important 
application of the power counting, one can also show that the interaction of 
the quark with ultrasoft gluons is suppressed by with respect to the non- 
relativistic limit described in the Schrodinger equation (|^). This means that 
up to NNLO0 we can simply neglect the effects of ultrasoft gluons — a great 
simplification because ultrasoft gluons cause Lamb-shift-type retardation cor- 

t The non-relativistic limit described in the Schrodinger equation {|^) is usually called 
"leading-order" (LO); corrections of order v" with respect to the non-relativistic limit are 
called N"LO. 



6 



rections that are quite difficult to calculate. For QQ productiop—they have 
never been fully calculated to date; only partial results exist so farE3. Another 
important issue is that the EFT gives an exact definition for the energy E in 
Eq. (|l|) , which is particularly important for the determination of heavy quark 
masses from experimental data. 



4 The Cross Section cr(e+e ^ QQ) Close to Threshold at NNLO 

To determine the total QQ production cross section in e+e^ annihilation in the 
non-relativistic region at NNLO in the EFT, we start from the corresponding 
expressions in full QCD {R = a{e~^e~ QQ)/a{e~^e~ — > 



■ Im 



i d^xe'''--{0\Tj^{x)j^iO)\0) 



■ Im 



-t{0\Tj,{q)r{-q)\0] 



(3) 



where eg is the heavy-quark electric charge and the cm. energy; j^(±q) 
are the electromagnetic currents that produce and annihilate the QQ pair with 



cm. energy 



. In the EFT these (external) currents are replaced by a sum 
of ^5*1 EFT currents with increasing dimension. For NNLO we need currents 
up to dimension 8 (i = 1, 2, 3): 

Mq) = ci - ^ Dfxjiq) + ■■■, 



U-q) = ci (x^<y^i){-q) - (x^^^{-h D)^){-q) + 



(4) 



where the Pauli spinors ip and x describe the heavy quark and antiquark fields 
in the EFT, respectively. The constants ci and C2 are Wilson coefficients that 
contain a perturbative series in originating from integrati|iiff|J3|Ut the off-shell 
d.o.f.'s. At NNLO ci has to be calculated at order al (Refs.llil't3) and contains 
only contributions from hard d.o.f.'s. For C2 the Born expression is sufficient. 
The dependence of ci on the regularization parameter is not written out. We 
note that only the spatial components of the currents contribute at the NNLO 
level. Using the expansions we arrive at the following expression for the 
cross section at NNLO 



KLoiq^)^^C^lra Ai{q^) 



"Q 



4 TT e?) r 
— -^C2lm A2{q^ 
6 rrin L 



(5) 



For simplicity we consider only the photon-mediated cross section. 
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where 



(6) 



A2{q^) = iz(0|(^Vx)(9) {x^a{-iD)^){-q) +h.c.|0) 



(7) 



are EFT current-current correlators and mq is the heavy quark pole mass. 
Using the EFT equations of motion for the heavy quark fields A2 can be 
related directly to Ai- 

A2^mQEAi. (8) 

The current-current correlators are also regularization-parameter-dependent 
and, in fact cancel the regularization-parameter dependence of Ci/2- Ai could 
be calculated from Feynman diagrams in the EFT. However, even in LO we 
would have to resum an infinite number of them (and we would have to carry 
out this rcsummation at NNLO). ft is more convenient to consider the four- 
point function (p = i\/q^, 0)) 



Gik,k',q') = (0 



p 







(9) 



that describes off-shell elastic scattering of a QQ pair in the EFT with cm. 
energy Integrating G over k and k' we arrive at Ai{q^) (up to a trivial 

factor) . Upon integration over ko and kg the equation of motion of G at NNLO 
is just the Schrodinger equation (|l|) modified by NNLO corrections.: 

■.2 



2mQ - 2(5toq) 



G{k,k';q^ 



'^'P V{k,p)G{p,k';q') 



(2^) 



(2 7r)3,5(3)(fc- fc')- 



(10) 



All terms displayed in Eq. (^^ are of order v"^ ~ al (LO) and higher. The 
potential V contains_the Coulomb potential shown in Eq. (|^), its one- and 
two- loop corrections 113, the order as/rriQ Breit-Fermi potential known from 
positronium, and an order o^./'itiq non-Abelian potential. This equation can 
be solved either perturbatively, starting from the known LO Coulomb solution, 
or exactly, using numerical methods. It is important that this calculation is 
carried out in exactly the same regularization scheme as was used to determine 
the constants C1/2. As mentioned before, we obtain Ai{q'^) by integrating 
G{k, k' ; q^) over k and k'. From the EFT we can also uniquely determine 
the exact definition of what is E in Eq. (^: in the pole mass scheme we find 
SruQ = in Eq. 
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5 The Heavy Quark Mass 



The pole mass definition seems to be tfie natural choice to formulate the non- 
relativistic effective theory that describes the QQ dynamics close to threshold. 
After all, the heavy quark pole mass is gauge-invariant and IR-finitelia. More- 
over, in the pole mass scheme the equation of motion for the non-relativistic 
QQ pair has the simple form of Eq. (|l|) (and also of (|l^) for Smq = 0), which 
is well known from the non-relativistic problems in QED. Intuition also seems 
to favour the pole mass definition, because close to threshold the heavy quarks 
only have a very small virtuality of order rriQv'^, i.e. they are very close to 
the mass shell. However, it is known that the use of the pole mass can lead 
to (artificially) |Laxge high order corrections, because of its strong sensitivity to 
small momentaEj. Quantitatively, the pole mass is ambiguous to an amount of 
order Aqcd . Thus the use of the pole mass is unfavourable if we intend to use 
the calculation for the QQ cross section at NNLO for a precise determination 
of the heavy quark mass. [At this point I emphasize that the notion of the 
"heavy quark mass" is not a physical one, because we can only define quark 
masses as parameters in the QCD Lagrangian, not as physical quantities.] In 
fact, the first theoretical calculations of the NNLO cross section for produc- 
tion close to thresholdEj'EZlEa were carried out in the pole mass scheme; they 
seemed to indicate a rather large instability of the location where the tt cross 
section at threshold rises, because this rise did not sepm to converge when 
comparing LO, NLO and NNLO calculations. (See Ref.Ej for a synopsis of all 
calculations that have been done for the total tt cross section at NNLO.) The 
situation was similar for the first determinations of the bottom quark mass 
from moments of the-tptal bb cross section at NNLO, where the uncertain- 
ties were quite large E2lEil if the pole mass was employed and the latter were 
estimated conservatively. 

The bottom line is that we should use a heavy quark mass definition that 
has an infrared sensitivity much smaller than the pole definition. In principle, 
the MS mass definition might be an ideal candidate because its definition does 
not depend on infrared physics. However, using the MS mass would mean that 
SniQ = itiq'^° —TTiQ is of order as and parametrically larger than any other term 
in the Schrodinger equation ( 10 ) . This means that the MS mass formally breaks 
the non-relativistic power counting. In order to find an adequate "threshold 
mass" definition, it helps to have a closer look at Eq. (|l^). Technically, the 
instability observed in the pole mass scheme originates from the potential V. It 
can be shown that the potential causes large corrections from momenta smaller 
than mqas at large orders of perturbation theory. This can be visualized by 
considering the small momentum contribution to the heavy quark potential 
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in configuration space representation for distances of the order of the inverse 
Bohr radius l/mQas'. 



Vc{q) exp{iqr) (11) 



(27r)3 

\q\<ti^Mta,s 

I (27r)3 

where Vc is the Coulomb potential in momentum space representation. At 



"^'^■14(9) + ..., (12) 



large orders of perturbation theory the RHS of Eq. (12) grows asymptotically 
like —fia^nl (Ref.EHl), where /i is the scale of as- On the other hand, it can be 
shown that the total static energy 2mQ°'° + Vc{r) does not contain these large 
corrections OEaS. Therefore, an adequate mass definition can be found by 
reabsorbing the r-independent small-momentum contribution of the Coulomb 
potential into the heavy quark mass. This procedure is notpUmnufi. Several 
different mass definitions have been suggestedin the literatureE2lc3L3 that fulfil 
this task. Let me mention only the IS massEZI in some detail. The IS mass is 
defined as one half of the mass of the perturbative contribution of a fictitious 
n = 1, ^Si QQ bound state, assuming that the quarks are stable: 

m'i~ml°'^ = l{l'S^\n\l'S,) = -^a^mg + ... , (13) 

where Ti. is the non-relativistic perturbative Hamiltonian. The NNLO expres- 
sion forjlJie RHS of Eq. (|l^), which is not displayed here, was first calculated 
in Ref. Ea. The IS mass is less sensitive to infrared momenta than the pole 
mass (its ambiguity is only of order Aqqj^/toq) and, at the same time, does 
not break the non-relativistic power counting. 

In order to relate the IS mass to the MS mass, some care has to be taken, 
because it is not a priori clear how to combine the corresponding perturbative 
series. The guiding principle is the proper cancellation of the asymptotically 
large perturbative corrections contained in the perturbative relations of the IS 
and the MS masses to the pole mass in each order of perturbation theory. It 
can be shown that the cancellation happensjproperly if a modified perturbative 
expansion, called the "upsilon expansion" o is used. In the upsilon expansion 
one has to consider terms of order a" in Eq. (|3|) as of order a"~^. Technically 
this unusual prescription comes from the fact that the coefficient of the order 
a" term in Eq. ([l3| ) contains the contribution X]"=o^ ^'^(^/'^s)/*'' which leads 
asymptotically to an additional factor l/a^. The upsilon expansion should 
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be applied whenever the IS mass is used outside the framework of the non- 
relativistic problems. As an examnle, it has been applied for the description 
of inclusive semileptonic B decays cB, where it leads to a very well behaved 
perturbative expansion. 

6 Top-Antitop Production Close to Threshold 

Top-antitop pair production at future lepton pair colliders is one of the most 
important applications of the EFT for non-relativistic QQ pairs discussed 
in this talk. The tt production close to threshold will provide a rich phe- 
nomenological environment where the top quark can be studied in quite a 
unique manner. Assuming the upcoming Linear Collider or a muon collider 
will spend around a 100 fb^^ at the top-antitop threshold, the top mass can 
be determined from the tt line shape with experimental uncertainties below 
50 MeV. The fact that the tt pair is produced in a colour singlet state is 
crucial to make this possible. Even the top quark width might be extracted 
from the line shape with an uncertainty of around 20%, assuming that the 
theoretical uncertainties are small. Assuming that the Higgs mass is around 
100 GeV also a fairly good determination of the top Yukawa coupling could 
be achieved. In addition, from studies of distributions one might get measure- 
ments of anomalous top couplings and further measurements of the top width. 
What makes the tt system so interesting theoretically is the fact that its total 
width (Tt w 1.5 GeV) is much larger than Aqcd- Thus, the top quark decays 
before it hadronizes, and the effects of non-perturbative effects should be small. 
One can show that the effective velocity of the top quarks close to threshold 
reads VeS = | + iTt)/mQ\ rather than y^E/niQ. This means that the hi- 
erarchy mt 3> mtVcS 2> mtv'^g > Aqcd holds even for cm. energies that are 
extremely close to the threshold and that non-perturbative effects should only 
represent relatively small corrections. 

In view of this prospect, it is obvious that a careful analysis and assessment 
of theoretical uncertainties in the prediction of the total cross section is manda- 
tory, in order to determine whether the theoretical precision can meet the ex- 
perimental one. Within the last two and a half years, considerable progress 
has been achieved in the calculation of NNLO corrections to the tt production 
process within the framework of the EFT for non-relativistic QQ pairs. For 
the total c.rop.p^i|±imJ-he complete set of NNLO QCD corrections have been 
calculatedlljOliJilj'Lilild. So far, for the total cross section the complete elec- 
troweak effects have only been implemented at NLO. The implementation up 
to NLO is quite easy: the only thing to do is to replace (y^ — 2 mg — 2 (5toq) 
by — 2mQ — 2 6mQ) + iTt in Eq. (|lo|). The width has to be counted 
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Figure 1: The total photon-induced tt cross section divided by the point cross sec- 
tion at the LC versus the cm. energy in the threshold regime at LO (dotted curves), 
NLO (dashed) and NNLO (solid) in the IS (left figure) and the pole (right figure) 
mass schemes for as{Mz) = 0.118 and fi = 15, 30, 60 GeV. The effects of beam- 
strahlung and initial-state radiation have not been included. The plots are taken from 
Ref.2''. 



as of order ~ v'^g and cannot be treated as a perturbation. In the left 
picture of Fig. |l| the result for the total tt production cross section is shown for 
ml^ — 175 GeV, where dotted, dashed and solid curves represent LO, NLO 
and NNLO results, respectively. As a comparison, the corresponding results 
are shown in the pole mass scheme in the right picture of Fig. |l|. It is conspic- 
uous that the normalization uncertainlifis are rather large, at the level of 20%. 
However, realistic simulation studies l3 have shown that these uncertainties 
do not severely affect the determination of a threshold mass such as the IS 
mass. A top mass determination with combined theoretical and experimental 
uncertainties of order or smaller than 200 MeV seems realistic. 



7 Other Applications and Conclusions 

There are a number of other applications of the EFT, which I will only discuss 
briefly. As mentioned in the introduction, moments of the total bb cross section 
in e~^e~ annihilation can be used to determine the bottom quark mass. These 
moments are defined as 

f dq^ ^e+e-^&6+x(g') 

"-y (<z)2("+i)cTe+,-^,+,-(<z2)- ^^^^ 

It can be shown that the effective velocity of the bottom quarks in the n-th 
moments is VcS = 1/y/n, which means that mf, > rrihVcff > TOfeU^ff > Aqcd 
holds as long as n is not chosen larger than about 10. Determinations of the 
bottom quark mass based on threshold mass definitions have been carried out 
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in Refs.c2l'OE3. The results are perfectly consistent and yield a value for the 
MS bottom mass mf,(mb) of 4.2 GeV, with an uncertainty of 60-80 MeV. 

It is clear that the non-relativistic EFT for heavy QQ pairs can also be ap- 
plied to two-body systems where the binding is caused by the electromagnetic 
forces. Apart from classic systems such as hydrogen, positronium or muonium, 
where the use of effective field theoretical methods has become increasingly 
popular within the last two years, hadronic atoms are a quite interesting field 
of application. Among the most interesting of these systems are the bound 
state of two charged pions (tt+tt"), called pionium, or muonic or kaonic hydro- 
gen. A number of experiments such as DIRAC at CERN or DEAR at DA$NE 
are already running or are under way, with the aim of measuring parameters 
of the chiral Lagrangian with unprecedented precision in order to provide tests 
of our understanding of chiral symmetry breaking. 

To conclude, one can say that the advent of EFT methods has revived the 
interest of many theorists in non-relativistic two-body systems. The strength 
of the EFT approach is that it represents a "theory" and not just a "method" . 
It explicitly uses the concepts of separation of scales, factorization and renor- 
malization, and it simplifies systematic calculations through a set of power 
counting rules that can be applied before any calculation is actually started. 
This has made non-relativistic heavy QQ systems accessible to a wider public. 
Today it is clear how to calculate NNNLO corrections to the total QQ cross 
section. This is a tremendous technical project, but it is clear conceptually. 
The calculation of these corrections might be important in obtaining a more re- 
alistic estimate of theoretical uncertainties in the non-relativistic expansion for 
and hh systems. An important conceptual development, which still has to be 
achieved, is the proper implementation of the instability of a heavy quark such 
as the top into the EFT framework. To study instability, the non-relativistic 
framework seems to be a much better setting than the full theory, because the 
non-relativistic EFT provides a natural separation of on- and off-shell heavy 
quark d.o.f.'s. 
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